Equipartition of rotational and translational energy in a dense granular gas 
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Experiments quantifying the rotational and translational motion of particles in a dense, driven, 2D granular 
gas floating on an air table reveal that kinetic energy is divided equally between the two translational and one 
rotational degrees of freedom. This equipartition persists when the particle properties, confining pressure, pack- 
ing density, or spatial ordering are changed. While the translational velocity distributions are the same for both 
large and small particles, the angular velocity distributions scale with the particle radius. The probability distri- 
butions of all particle velocities have approximately exponential tails. Additionally, we find that the system can 
be described with a granular Boyle's Law with a van der Waals-like equation of state. These results demonstrate 
ways in which conventional statistical mechanics can unexpectedly apply to non-equilibrium systems. 



PACS numbers: 45.70.-n, 45.50.-j, 45.70.Mg 



Many processes operate far from equilibrium and there- 
fore cannot be described within the framework of classical 
thermodynamics and statistical mechanics. One of the sim- 
plest non-equilibrium systems is the granular gas, composed 
of driven, colliding macroscopic particles [1 1. Such a system 
differs from an ideal gas in two important ways: energy in- 
put is required to maintain a gaseous state, and steric effects 
arise due to finite particle size. Granular gases can exist in 
a non-equilibrium steady state (NESS) in which energy flows 
into the system from an injection mechanism and out of the 
system via dissipation. Studying granular gases allows us to 
examine the extent to which these systems can be described 
using the techniques of conventional statistical mechanics. 

While there has been a great deal of success describing di- 
lute granular gases using a Chapman-Enskog expansion of the 
Boltzmann equation pi ^31, deviations from ordinary gases 
arise as the density of the gas increases iH. For example, 
granular gases exhibit non-Gaussian velocity distributions |4 - 
|8J, non-equally partitioned kinetic energy (temperature) in 
mixtures of different particles iH |9l (TO), and a breakdown 
of molecular chaos due to velocity correlations (TT\. These 
behaviors, together with simulations in which strong coiTe- 
lations between translational and rotational velocities are ob- 
served [121, raise the question of how such systems partition 
energy (or momentum) between rotational and translational 
modes. In this Letter, we present an experiment which, re- 
markably, exhibits equipartition between the rotational and 
translational energies. 

To connect macroscopic state variables (pressure, volume, 
energy) with microscopic dynamics (translational and rota- 
tional particle velocities), we perform experiments on a dense 
granular gas composed of a single layer of bidisperse disks 
lfT3l [T4ll . The particles float on an air hockey table and are 
driven by impulses from bumpers at the boundaries, produc- 
ing a fluctuating but stationary state. In contrast with the many 
previous experiments [4- -, 11 J which inject energy by global 
agitation, we drive the system via particle-scale bumpers at 



the boundaries. We find that the distributions of the particle 
velocities are exponential in character and that they are re- 
markably unaffected by modifications to the packing density, 
particle interactions (friction/restitution) or the boundary con- 
dition (constant pressure or constant volume). Because the 
time-averaged energy of the particles is constant, we are able 
to identify a gas-like equation of state (EOS) with a van der 
Waals correction, similar to prior observations in a monodis- 
perse system of vertically-vibrated spheres 1151 . This EOS 
points to the utility of thermodynamic-like approaches, even 
for far-from-equilibrium systems. 

Experiment — We perform experiments on bidisperse 
disks (petii dishes) which float on a horizontal air table (see 
Fig.[T^; further description in lHUHH). We use a high speed 
camera to track the translational and rotational trajectories at a 
time-resolution sufficient to resolve individual collisions. The 
system is driven on three of four sides by impulses from elec- 
tromagnetic bumpers; random bumpers are triggered every 
0.1 s. As shown in Fig.[T]3, the injection of kinetic energy is an 
iiTegular process, with intermittent strong heating followed by 
approximately exponential cooling. This irregularity is likely 
due to transient force chains lfT6l [17 1 which happen to termi- 
nate at a firing bumper and can therefore efficiently transmit 
impulses. The system is constrained on the fourth side by a 
piston (mass = 439 g) which also floats on the air table. The 
position of the piston can be fixed to create a constant vol- 
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FIG. 1: Color online, (a) Schematic of experiment, top view, with 
bumpers at sides. The piston provide CP boundary conditions or be 
fixed in place to create CV. (b,c) Timeseries of average kinetic energy 
per particle at two values of (CV boundary conditions). 
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ume system (CV), while a constant pressure system (CP) is 
produced by suspending masses on lines which run from the 
piston to pulleys on the opposite side of the table. The two- 
dimensional pressure exerted on the system ranges from 0.01 
N/m < P < 0.1 N/m. We primarily investigate amorphous 
packings with packing density 0.61 < <j) < 0.79 (147 to 186 
particles). These (f> lie below the point at which a jamming 
transition would occur in a static system |18, 19 1. The up- 
per limit of the accessible (j> range is set by the pressure at 
which out-of-plane buckling occurs, while the lower limit is 
set by the requirement that (p must be high enough for firing 
bumpers to contact particles. In addition to amorphous pack- 
ings, we have also examined crystalline packings in which the 
system is prepared with all of the small particles on the +y 
side of the table and all of the large particles on the —y side. 

We prepare the system with N = Ns + Nl particles so 
that the area covered by large and small particles is equal. 
Large particles have radius = 41.8 mm and small par- 
ticles rs — 27.9 mm; since (rL/rs)'^ = 2, this requires 
Ns = 2A^i. For bare petri dish contacts, the coefficient 
of friction is /i = 0.5 and the coefficient of restitution is 
e = 0.33 ± 0.03 ||T3l . The friction and restitution coeffi- 
cients can be adjusted by stretching an elastic band around 
each of the particles so that the diameter remains unchanged; 
this results in // = 0.85 and e' = 0.51 ± 0.07 HU. (All 
variables describing the banded particles are denoted with the 
prime.) The mass m and moment of inertia / for bare dishes 
are rriL = 8.13 ± 0.20 g and II = 1.3 ± 0.1 x lO"* g 
mm^ and ms = 3.52 ± 0.12 g and Is = 2.5 ± 0.2 x 10^ 
g mm^. With the addition of rubber bands, m increases by 
0.44 g, and / increases to Ig ~ 2.8 ± 0.2 x 10^ g mm^ and 
= 1.4 ± 0.1 X 10* g mm^. The moment of inertia is mea- 
sured by rolling particles down an inclined plane. 

Each particle is tagged with a unique nine-dot pattern which 
allows the identity, position and angular orientation to be de- 
termined; in addition, we track the location of the piston to 
measure the system area A. We film the experiment using 
Phantom V5.2 high speed camera and measure the transla- 
tional {vx, Vy) and rotational (w) velocities from the trajectory 
of each particle. A key benefit of unique labels is to eliminate 
error caused by confusing the identities of two particles |20|. 
Images are taken at a frame rate of 40 Hz, and a 5-point box- 
car average is applied along the trajectories to reduce noise. 
Sampling at a lower rate causes the velocity distributions to 
narrow, while taking images at a higher rate leaves the ve- 
locity distributions unchanged. Measurements of stationary 
particles indicate that the v resolution is ±2 mm/s, while the 
Lj resolution is ±0.4 rad/s. The 40 Hz frame rate allows each 
measurement run to last 100 s (4000 frames); 1 — 4 such runs 
are collected for each choice of (p, P, the boundary condition 
(CP, CV) and the particle properties (bare, banded). We find 
that stationarity of statistical distributions is achieved within 
several seconds. 

Velocity distributions — As shown in Fig.|2^,b, the prob- 
ability distributions of the translational velocities Vx.y are ap- 
proximately exponential in character, and have the same mean 
speed \ vx^y \ ~ 20 mm/s. Several observations suggest that the 
broader-than-exponential tails of 'P{v) are due to high-speed 
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FIG. 2: Color online. Translational velocity distributions (a) Vivy), 
with {\vy\) = 19.9 mm/s (b) V{v:c), with {\vx\) = 19.4 mm/s and 
(c) 'P{vy) for particles which are at least 3r_t from a boundary, with 
{\vy\) = 19.3 mm/s. (d) Rescaled rotational velocity distribution 
Vlrui), with {r\uj\) = 20.9 mm/s. Colors indicate conditions of run 
(large, bare, amorphous, unless otherwise specified.) 

particles near the bumpers. First, V{vx) is asymmetric: the 
(— x) side of the distribution is narrower because there are no 
bumpers along the piston. Second, if we restrict 7'(tij^) to only 
those particles more than ir^ from the boundaries, a more 
exponential form is recovered, shown in Fig. [2]:. The high- 
velocity tails for the small particles are slightly wider than for 
the large particles, which suggests that the small particles ac- 
quire a slightly larger velocity from the bumpers. 

The angular velocity distributions V{io) are also approxi- 
mately exponential. In addition, the V{'jj) distributions col- 
lapse to a single curve when rescaled by the particle radius 
(see Fig. |2]i), suggesting that the system prefers rolling con- 
tacts. Measurements of the relative velocity at each contact 
confirm this preference. The shapes of all three Vivx, Vy,Ld) 
distributions are very robust: switching from CP to CV bound- 
ary conditions, varying cj), changing ^ and e, and ordering the 
system have little effect on any of the velocity distributions. 
In addition, the mean values of (l^a;!), (|t;j,|), and {r\uj\) are 
all similar. As we will examine below, this similarity results 
from both equipaitition and the way mass is distributed within 
the particles. 

Equipartition — The time average of total kinetic energy 
of the system, E = {E{t)), is the sum of the translational 
(Eti) and rotational (i?iot) energies of all the particles (i — 
1..N) in the system: 

^ = E l^Mcc) + E I'^My) + E l^^i^')- (1) 

i i i 

The rms velocities, (•), are averaged over time. In equilib- 
rium thermal systems, equipartition requires that each degree 
of freedom (DOF) have the same energy; however, no such 
requirement exists for non-equilibrium systems. As shown in 
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FIG. 3: Color online, (a) Time-averaged ratio of the translational ki- 
netic energy to the total kinetic energy (Etr/E) averaged over 100 s. 
Dashed line is 2/3, corresponding to equipartition. (b) Average ki- 
netic energy per particle as a function of (f). Error bars are typical 
standard deviation. Each data point indicates the average for 25 s of 
data (CP) or 100 s (CV). Color of symbols indicates P (see Fig.|4] 
for key). 



Fig. [2] we observe that the large and small particles have the 
same indicating that translational energy is not parti- 

tioned equally between the two species. However, the sys- 
tem does partition the energy equally between the two trans- 
lational and one rotational DOF: as shown in Fig.|3^, the value 
of (Etr/E) is approximately 2/3. This observation holds for 
both bare and banded particles and for both CP or CV bound- 
ary conditions, although the ratio Eti/E decreases slightly 
with increasing <j) (or P). 

In Fig. [2jl we observed that (r^lojil) = {rs\u>s\)- How- 
ever, many different pairs of V{ijJs) and Viuj^) can satisfy 
this requirement. An additional constraint is provided by the 
equipartition of Etr and iJrot : 



(Nsms + NLmL){v, 



x,y/ 



NrL 



(2) 



The moment of inertia can be written as li — airriirf, 
where a describes how mass is distributed over the particles: 
as = 0.97, and = 0.94. Eq. [2]can be satisfied for any 
two values of aL and as when (r^lwil) = {rs\ijJs\) and 
T^iys) = 7^(i'l). The observation that Et^/E is systemati- 
cally lower for the large particles follows from the observa- 
tion that a^ < as- Because the two have the form 
g-a|a;i|/ri ^^^^j ^-a\u:s\/rs ^ j^is pair of distributions can fulfill 
the geometrical preference for rolling contacts for any value 
of a selected by the system. The system is observed to choose 
a such that there is equipartition of rotational and translational 
energy at each 4>. 

Equation of state — In equilibrium statistical mechanics, 
the total kinetic energy (Eq. [T]) corresponds to the tempera- 
ture: E = '^ksT. In our expeiiment, we observe that even 
though E{t) fluctuates significantly (see Fig.[lJ)), the average 
per-particle value E/N = {E{t))/N is a constant indepen- 
dent of (p, material, boundary condition, and degree of spatial 
ordering (see Fig.jSj)). This value is not universal, but depends 
on such details as bumper firing frequency, piston inertia, and 



number of particles. In analogy with a thermal bath at con- 
stant T, the bumpers at the boundary create a granular gas 
with constant E. It is therefore natural to inquire whether a 
granular version of Boyle's Law {PV = const.) holds. 

In two dimensions, the pressure exerted by the piston is 
given by a force per unit length, so that P = F/L and PA has 
units of energy. As has been done for experiments on vibrated 
granular crystals |15|, we apply a van der Waals-like correc- 
tion to account for the excluded area. The free area available 
to the particles is the area under the piston A{t), less the area 
occupied by the particles and the interstices too small to ac- 
commodate another particle. This excluded area is given by 
A* = {NsTTTg + 7Vi7rr|)/ (p*, where (/)* is a maximum pack- 
ing density. In an ideal gas PA = E, but in dense systems the 
compressibility increases due to the relative inaccessibility of 
nearby configurational states. These steric effects introduce 
an additional factor (1 + x) lED> where the compressibility 
factor X is proportional to 1 + e and to the pair correlation 
function g{r) evaluated at the paiticle radius. Including all 
three considerations, we consider the following EOS: 



PiA-A*)^E{l + x). 



(3) 



In spite of the large temporal fluctuations of E{t), E — 
{E{t)) is constant, as shown in Fig.jsj?. We test Eq.|3]by vary- 
ing P and observing A{t) for a fixed number of either bare 
or banded particles. We find that there is indeed a linear rela- 
tionship between A and P (Fig. |4]i, where the intercept corre- 
sponds to the excluded area A* . For P < 50 mN/m, the fluc- 
tuations in the area under the piston are large and the system 
does not behave as though it is in a steady state on the time 
scale of the experiment. From a least-squares fit weighting 
each point equally, we determine 0* = 0.79 ± 0.02 and * = 
0.79 ± 0.04 for the bare and banded particles respectively. For 
comparison, previously measured 1 13] values of static random 
loose packing estimate 0rlp = 0.81 for both bare and banded 
particles. For the crystalline system, 0*i.yst = 0.81 ± 0.02, 
which is also less dense than the hexagonally close packed 
value (/)HCP = 0.91. Therefore, the measured (j)* do not cor- 
respond to the van der Waals-like maximally dense static state 
(this would corresond to either random-close-pacing (j)RCP or 
4>HCp)^ but to some kinetically-determined value. 

The slope of the equations of state plotted in Fig. |4] cor- 
respond to 1 + x; their deviation from a value of 1 corre- 
sponds a measurement of the compressibility. For bare parti- 
cles, 1+x = 1.3±0.1, for banded particles 1+x' = 3.2±0.8, 
and for crystalline bare particles 1 + Xcryst = 0.7 ± 0.1. 
Although X is predicted to increase as the system becomes 
denser [21 1, we instead find that x is approximately constant 
for all the values of P that we examined. This result is un- 
expected, especially given that experiments on a single layer 
of vertically-shake spheres ifTSll indicate that x rises with in- 
creasing P as expected. 

Conclusions & Discussion — We have tracked the transla- 
tional and rotational motion of a dense granular gas of bidis- 
perse particles driven by bumpers at the boundary of a hori- 
zontal air table. The states generated have constant P, A, and 
E, which motivates treating the system within the framework 
of statistical physics. The velocity distributions are approxi- 
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FIG. 4: Color online. Data from CP experiments, corresponding to 
Eq.[3] Symbol color indicates P; symbol shapes as shown in Fig. [3] 

mately exponential and are independent of the packing den- 
sity, particle properties, and degree of spatial ordering. The 
exponential character of these distributions is consistent with 
both horizontal quasi-2D free cooling experiments [4J and 
quasi- 2D vibration experiments in microgravity fS] and ||8]. 
Our system appears to experimentally realize a state in which 
rare, high-energy injections produce a well-mixed steady state 
ll22l . although simulations of this driving mechanism find 
power law tails in the velocity distributions, while we ob- 
served nearly-exponential tails. Although energy is not par- 
titioned equally between the large and small particles, we do 
observe equipartition between the two translational and one 
rotational DOR Finally, the system exhibits a van der Waals- 
like equation of state. These findings indicate that while gran- 
ular gases can exhibit some properties of equilibrium statis- 
tical mechanics (equipartition of rotational and translational 
energy, Boyle's Law), they do not necessarily do so (the veloc- 
ity distributions are not Gaussian and energy is not partitioned 
equally between the large and small particles). The observa- 
tion of a thermodynamic -like EOS for such a small system 
provides encouragement that similar descriptions may be pos- 
sible in the many boundary-drive granular systems where only 
a mesoscopic portion of the system, typically < 10 particle di- 



ameters within shear bands, is set into motion. 

We find that the translational velocity distributions are the 
same for both the large and small particles; consequently, 
energy is not partitioned equally between the two particle 
species. Other experiments on ensembles of particles made 
from different materials |4 , 9 , 10 1 indicate that energy is not 
partitioned equally between particles with different properties. 
The absence of equipartition amongst particles with different 
restitution coefficients has been explained within the context 
of Chapman-Enskog expansions [24 1. In addition, simulations 
Il25ll26l predict that the way energy is distributed between ro- 
tational and translational modes depends on /i and e; however, 
for the values used here, equipartition of Eti and i^rot is not 
predicted. 

While we observe equipartition between the rotational and 
translational degrees of freedom, it is notable that no such 
equipartition is observed in dilute granular gases with fully 3D 
motions |27|; this raises the question of whether the density 
or the dimensionality is the key difference which promotes 
equipartition. Additional constraints arise due to density: par- 
ticles most easily translate in concert with their neighbors and 
rotational motion is governed by interparticle friction. That 
V{vx^y) are the same for both the large and the small particles 
in our 2D system indicates that translational motion requires 
cooperation. In addition, that 'P{vrUj) are the same for both 
the large and small particles indicates a preference for rolling 
contacts. The latter constraint does not apply in 3D because 
translational motion can occur along the axes of rotational mo- 
tion. Once the preference for rolling contacts is fulfilled, our 
2D system still must select how to distribute kinetic energy 
amongst the rotational and translational degrees of freedom. 
Unexpectedly, our findings suggest that equipartition, while 
not guaranteed for a non-equilibrium system, is used to select 
how the rotational and translational energy is distributed. 
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